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Abstract 

In this paper, we study finite dimensional approximations of Kolm-Sfiam models, which 
are widely used in electronic structure calculations. We prove the convergence of the finite 
dimensional approximations and derive the a priori error estimates for ground state energies 
and solutions. We also provide numerical simulations for several molecular systems that support 
our theory. 
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1 Introduction 

Density functional theory (DFT) is a theory of many-body systems and has become a primary tool 
for electronic structure calculations in atoms, molecules, and condensed matter [TBI [HI [HI [Ml [25l 
I26j . The most widely used is the Kohn-Sham model, in which a many-body problem of interacting 
electrons in a static external potential is reduced to a tractable problem of non-interacting electrons 
moving in an effective potential. The purpose of this paper is to analyze the finite dimensional 
approximations of Kohn-Sham models so as to provide a mathematical justification for both the 
directly minimizing energy functional method |241 \T7\ and the variational optimization method (i.e. 
solving the Kohn-Sham equation self-consistently) [23] and some understanding of several existing 
approximate methods in modern electronic structure calculations. 

Throughout this paper, we restrict our mathematical analysis and numerical simulations to non- 
relativistic, spin-unpolarized models. In the pseudopotential setting, the ground state solutions of 
the Kohn-Sham model for a molecular system can be obtained by minimizing the Kohn-Sham 
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energy functional 



E{{(|>^}) = o V / \yMx)\^dx+ Vioc{x)p{x)dx + S2 M^)VniM^)d: 




\x 




(1.1) 



with respect to wavefunctions under the orthogonahty constraints 




where A'' is the number of valence electrons in the system, p = is the electron density, 

Yioc and Vni are the local and nonlocal pseudopotential operators respectively, that treat the core 
electrons and the nuclei as a unit and represent the interactions on the valence electrons [23], and 
f (p) denotes the exchange-correlation energy per unit volume in an electron gas with density p. The 
Euler-Lagrange equation corresponding to this minimization problem is the so-called Kohn-Sham 
equation: find \i e K, e H^i^^) = 1, 2, • ■ • , iV) such that 



where is the effective potential relative to the last four terms in energy functional (|1.1[) . 

This is a nonlinear integro-differential eigenvalue problem, and (|1.2p is often called self-consistent 
field (SCF) equation as to emphasize the nonlinear feature encoded in V'e//({(/>i}). It is assumed 
in most of the simulations that the ground state solutions can be found by occupying the lowest 
eigenstates of Kohn-Sham equation (jl.2p . It is not known whether the assumption is true, but it 
seems to be most often the case in practice. 

The main difhculties of numerical analysis for Kohn-Sham models lie in what we have to either 
handle the global minimization problems whose energy functionals may be nonconvex or deal with 
the nonlinear eigenvalue problems whose eigenvalues may not be nondegenerate. To our best 
knowledge, except for the very recent works of Cances, Chakir, and Maday [B] and Suryanarayana 
et al [29], there is no any other numerical analysis for Kohn-Sham models in the literature. We 
see that the numerical analysis of Kohn-Sham models is crucial to understand the efficiency of the 
numerical methods widely used in electronic structure calculations. Under a coercivity assumption 
of the so-called second order optimality condition, [6] provided numerical analysis of plane wave 
approximations and showed that every ground state solution can be approximated by plane wave 
solutions, and [29] gave the convergence of ground state energy approximations based on finite 
element discretizations only. In this paper, we shall present a systematic analysis for a general finite 
dimensional discretization and prove that all the limit points of finite dimensional approximations 
are ground state solutions of the system, and every ground state solution can be approximated by 
finite dimensional solutions if the associated local isomorphism condition is satisfied. We provide 
not only convergence of ground state energy approximations but also convergence rates of both 
eigenvalue and eigenfunction approximations. We point out that the local isomorphism condition 
should be very mild and is indeed satisfied if the second order optimality condition is provided. 

Besides the Kohn-Sham models, there is another approach in DFT that is not so popular and 
is called of orbital- free DFT [lOl [31] , in which approximate functionals in terms of electron density 
alone are used for the kinetic energy of the non-interacting system and only the lowest eigenvalue 
needs to be computed. There are several related works on its convergence analysis [H [191 ISS] [33] 
and a priori error estimates [5l [6l [9] . 





(1.2) 
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This paper is organized as follows. In the coining section, we give a brief overview of the 
Kohn-Sham models and some preparations. In Section [31 we derive the existence of a unique local 
discrete solution under some reasonable assumptions. In Section 21 we prove the convergence of 
finite dimensional approximations of the ground state solutions with quite weak assumptions and 
derive the error estimates of ground state energy, ground state eigenfunctions and eigenvalues. 
In Section |S1 we present some numerical results that support our theory. Finally, we give some 
concluding remarks. 



2 Preliminaries 

Physically, the Kohn-Sham model is set over R'^. But in a lot of computations, M.^ may be replaced 
by some polyhedral bounded domain J7 C M'^, for example, a supercell for crystal or a large enough 
cuboid for finite system, which is reasonable since the solution of (II. 2p always decays exponentially 
[H [TSl I2H]- Thus we study numerical analysis of finite dimensional approximations of Kohn-Sham 
equation as follows: 

(-iA + K//({0J))0» = A,0, in^!, i^l,2,---,N, 

(2.1) 

n 

with the Dirichlet boundary condition (f>i — on dfl for finite systems and periodic boundary 
conditions for crystals, where C M"^ is a polyhedral bounded domain. 

We shall use the standard notation for Sobolev spaces VF^'^(Sl) and their associated norms and 
seminorms, see, e.g., [IT]- For p = 2, we denote H%Q) = W^^fl) and i?o(^) = {v & H^{n) : 
v\dn = 0}, where v\gn = is understood in the sense of trace, |1 ■ ||s,n = || • IU,2,n, and (■, ■) is the 
standard inner product. The space Y*, the dual of the Banach space Y, will also be used. For 
convenience, the symbol < will be used in this paper. The notation A < B means that A < CB 
for some constant C that is independent of the mesh parameters. 

Given ci £ R and p, C2 S [0, oo), we define 

^(p, (ci,C2)) = {/: 3 ai,a2 £ K such that cit^' + ai < J {t) < 02^^ + a-i V i > O}. 
For K £ R^^^, we denote its Frobenius norm by |k|. We consider the functional spac^l 
'H = {Hl{n)f ^{{c^^,<j>2,--- c^,^Hl{n) (j = l,2,.-. ,iV)}, 

/ N \ 

which is a Hilbert space associated with the induced norm ||$||i,o ~ j ^dl'/'illo.n + ll^'/'j llo.o) ) 

N 

and inner product (V*, V*) = ^(Vc?!),, V-^j) for $ = (^i, 02, • • • ,<j)N),^ = (■0i,V'2,- • ' , V'Jv) e U. 

i=l 

For simplicity of notation, we will sometimes abuse the notation by 

/ N \ . ^ > l/p 



^In fact, our theory also applies to space "H = {H^{Q))^ , where Q is the unit cell of a periodic lattice TZ of ] 
and H^{n) = {v\q : v e Hj^^iW^) and visTL- periodic}. 
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for subdomain uj C ft and $ = (0i,02,-'- 74>n) G For any $ = ((/)i,02,--' t4>n),'^ = 



N 



(V'l, • ■ • J V'Jv) G we define p$ = ^ l^ij^ and 



1=1 

AT 



\Jn J i,j=i 

In our discussion, we shall also use the following three spaces: 

= {M e R^""^ : = M}, A^^""^ = {M e R^><^ : = -M}, 

and 

Q = e : = /^^^}. 
We may decompose "H as a direct sum of three subspaces [HI [35] : 

H = 5$ ® ^$ ® 7i 

for any $ G Q, where 5ce, = $5^^^, = and %, ^ {"^ e H : '^^<S> = G M^><^} . 

2.1 Kohn-Sham models 

In the most commonly setting of local density approximation (LDA) [23] , the associated Kohn-Sham 
energy functional of (|2.ip is expressed as 

f ( ^ I ^ \ I 

•^^ \i=i i=i / 

for $ = 02, • • • ,'/'Af) G H, where VJoc is a smooth local pseudopotential, y„; is the nonlocal 
pseudopotential operator (see, e.g., [23]) given by 

M 

with Qj G L'^{Vl){j = 1, 2, • • • , Af), D{p<i,, p$) denotes electron-electron coulomb energy with 



nJn \x~y\ 

and £{t) is some real function over [0,cx)). We may assume that V/oc G i^(ri). We see that the 
function £ : [0, oo) — > M does not have a simple analytical expression. In applications, we shall use 
some approximations to f , for which we shall make the assumption that £{t) G .^(3, (ci, C2)) with 
Ci > or £[t) G ,^(4/3, (ci, C2)) that is satisfied by most of the approximations. 
First of all, we have 



Proposition 2.1. Functional (12.21) is invariant with respect to unitary transformations, i.e. 

E{^) = V $ G Q 

for any matrix U — {uij)fj^i G O^^^ , where O^^^ is the set of orthogonal matrices. 
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Using similar arguments in [5], we obtain that is bounded below over Q. More precisely, 

we have 

Proposition 2.2. There exist constants C > and 5 > such that 

E{^)>C-^\\^\\l^^^-h VJ-gQ. (2.3) 

To prove the convergence of the numerical approximations, we need the lower semi-continuity 
of the energy functional in the weak topology of whose proof can be referred to [S]. 

Proposition 2.3. // converge weakly to 5' in "H, then 

< liminf 

The ground state energy of the system is the global minimum of E{'^) in the admissible class 
Q and we shall study the following minimization problem 

inf {£;($):$ G Q}. (2.4) 

The existence of a minimizer of (|2.4p can be found in [31 [201 US or by similar arguments to that 
in the proof of Theorem 14.11 We see from Proposition 12.11 that if $ is a minimizer of (|2.4p . then 
G Q is also a minimizer for any U G O^^^ . Note that the uniqueness of a minimizer of (12.41) is 
open even up to an orthogonal transform since the energy functional may not be convex for almost 
all systems of practical interest. Therefore, we need to define the set of ground state solutions as 
follows 

= I $ G Q : E(^) = min E(^) 

We see that a minimizer $ — (0i, 02, • • • t4>n) of (12. 4p satisfies the associated Euler-Lagrange 
equation: 

N 

{A^M = (^Ay</)j,«) ^veHliSl), f = l,2,-.- ,7V, 

3=1 (2.5) 

n 

where is the Kohn-Sham Hamiltonian operator given by 

A^ = -iA + Vioc + Vni+ f f^dy + £'(p$) (2.6) 
2 Jn\- -y\ 

with the Lagrange multiplier 

A-(A.,)f,.^i = f / 0,^*0.) . (2.7) 

We define the set of ground state eigenpairs by 

e = {(A,$) G R^""^ X Q : $ G e and (A, $) solves ([23])} . 

Proposition 12.21 and (|2.7p imply that the ground state solutions are uniformly bounded 

sup (||$||i.n + |A|) <C (2.8) 
(A,*)Ge 
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for some constant C. 

To obtain the a priori error estimates of the finite dimensional approximations, we shah represent 
Kohn-Sham equation in another setting. Define 

with the associated norm ||(A, = |A| + ||$||i,o for each (A, <J>) £ Y. We may rewrite (|2.5I) as 

a nonhnear problem as follows: 

i^((A,$))=OGr*, (2.9) 

where F : Y ^ Y* is given by 

N N N „ 

(F((A, $)), {^,r)) = (^^-^^ - E ^u-^.'T.) + E / - (2.f0) 
j=i j=i 

with r = (71,72, • • • ,7w) G H and X = (x^,)5=i S K^^""^- 

The Frechet derivative i^,-'^ of at (A, ^) -.Y ^Y* \s defined as 

N AT „ 

= {C'^{K<^)^,T) - J2 + E / {i^^^J + M^J) V (p,*),(x,r) G r,(2.ff) 

where 

TV 



(/::,(A,$)*,r) = ii?"($)(vi/,r) - ^ (a,,Vj,70 



^2 

TV Af JV 

-{Y,Kj^,,li)+{2<j,,£''{p^)Y,<j,,ij,,^i)+Y.'^D{(l),^,,(l),-li)) (2.f2) 

for ^ = (^^1,^2, • • • , Vw) e and = (Ai„)f,=i 6 K^'^^. 
2.2 Basic assumptions 

The analysis of finite dimensional approximations will be carried out under certain assumptions, 
which are stated as follows 

Al \£'{t)\ + \t£"{t)\ e ^(pi,(ci,C2)) for some pi G [0,2]. 

A2 There exists a constant a G (0, 1] such that \£"{t) \ + \t£"'{t)\ < 1 + V i > 0. 

A3 If (A, <i>) is a solution of (12. 5p . then £'$(A, <i>) is an isomorphism from 7$ to 7$, namely, 
there exists a positive constant 7 depending on (A, $) such that 

mt sup -f— Ti rjT^Ti > 7. (2-13) 

*er*re7i ||«'l|i,o||r||i,n 
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We see that Assumption A2 implies Assumption Al and the commonly used and LDA 
exchange-correction energy satisfy Assumption A2 [F, "S; . We shall mention that the above assump- 
tions are necessary for the a priori error estimate, but none of them will be used in our convergence 
analysis of finite dimensional approximations (in Section HTTj) . 

Remark 2.1. It is open whether Assumption A3 holds for all Kohn-Sham models, though it may 
hold for semiconductors and "closed shell" atoms and molecules. We see that the following assump- 
tion 

{C'^{K<^)-^,-^)>^\mln V*er$, (2.14) 

which implies (|2.13[) . is employed in \ 27}j . Note that {2.1^^ is equivalent to (j2.13p when (A, <i>) is 
the ground state solution of (j2.5p . 

The following lemma will be used in our analysis of the local uniqueness of discrete solution. 

Lemma 2.1. Let ?/i = (Ai, $i) and y2 ^ {■^^2^ '^2) satisfy \\yi\\Y + \\v2\\y < C ■ If Assumption 
Al is satisfied, then there exists a constant Cp depending on C such that 

\\F{vi)-F[y2)\\<CF\\vi-y2\\Y ^VuV2eY. (2.15) 
Moreover, if Assumption A2 is satisfied, then there is a constant C'p such that 

\\F;,,-F'yJ<C'p{\\vi-y2\\^ + \\vi~y2\\l) ^yi.V2eY. (2.I6) 

Proof. To prove (|2.15p , it is sufficient to show that 

(^$,$i-A$,$2,r) <c||$i-$2lli,n||r|li,n vtgh, (2.17) 

which together with (|2.10p indeed implies (|2.15p . Using the Holder inequality and the Sobolev 
inequality, we have for i = 1,2, ■ ■ ■ , N that 

((-i A + Vloc)(t)l,^ - (-i A + Vioc)02,«, «) 

^ 2"'^^'* ~ '/'2,i||l,n||w||l,0 + ||V;oc||o,n||01,i - 02,i||o,3,n||t'||o,6,n 
< Ui,. - (l)2A\i,n\\v\\i,n yveH^ifl) 

and hence 

((-iA + v^<„)<fi-(-iA + ^<„)<f2,r) < l|<j>i-<j>2l|i,o||r||i,o vreH. (2.I8) 

Due to 

N M 
i=l j=l 

we obtain 

JV 

iVni^i~Vni^2,T) <Y,Uh^-^2,^\WnH\Wn < |1 $1 - $2 1 1 i,o || F 1 1 i^o V T G (2.19) 
Obviously 
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when pi = in Assumption Al. If Assumption Al is satisfied for pi e (0,2], then there exists 
G [0,1] such that 

N 

(£'(p*i)*i-^'(p*J'i>2,r) = V / (£'(p$J0i,,-f'(p*J</)2,.)7» 

N 



N 

< X! + '^&£"iPi)\\o,3/pi,n\\(t>l,i - 02,i||o,6,f2||7i|lo,6/(5-2pi),a 

1=1 

N 



< 



E ll^«lla3,oll</'i,. - hdiMkn < 11*1 - 'i>2||i,o||r|li,o, (2.20) 



1=1 



where ^ = (^1,^2, • • • An) with — 5i4>i^i + (1 — 5i)4>2,i, and the Holder inequality, the Sobolev 
inequality, and the fact 

llpdlo.3.o < ||p$Jlo,3,o + I1p*J1o,3,o < ll1>i|l?,o + \\'^2\\ln < 

are used. 

For Coulomb term, we obtain from the Young's inequality and the Holder inequality that 
\\r-^ * {p^, -p<i>J||o,oo,o < ||?-"^|lo,allP<i>i -P'S>2\\o,n < ||r"io^s^||$i - $2||i,n, 
where Q — {x ^ y : x,y (z Q}. Since 

n 

-1 , „ J, ^„, I / „-i 



(r * p^Ji(pi,i - (t>2,i)v + / r * (p*i - p$J(?!)2,iW 
/f2 Jo 

< \\r^^ * ||o,oo,o||<?!'i,j - 02,j||o,o|k||o,o + \\r^^ * (p*i - P$2)llo,oo,o||02,»||o,o||w||o,o 

< ||</'l,z-02,^||l,oM|l,O+||*l-*2||l,o|k||l,O ^veH^ifl) 

holds for i 1, 2, • • • , A^, we have 

((r-i*p$J$i-(r-i*p$,)$2,r)< ||$i-$2||i,o|lr|li,f2 VTG-W. (2.21) 

Taking (|27T8)) . ((219)) . ((2?20| . (|23T|) and definition into account, we then arrive at (l2T7l) . 

If Assumption A2 holds, then following [6] Lemma 4.5] we obtain for ^I^ = (■01,-02, ■ ■ • , ^w), T = 
(7ii72, • • • ,7Af) e ^ that 

. .1 iV AT 

|(f'(p$jvi/,r)-(£'(p*j«',r)| - / / 2£"{p^^,^){J2Mt){c|^l.^-h.i))C£'f''^^^^d^ 

•^f^ -^o »=1 z=l 

(2.22) 
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and 



N 



N 



N 



N 



O Jo 



j=i i=i i=i 

N N N 

1=1 1=1 

E - 02,.)) ( E ( E '^■'(^)tO 



< 



n j_ "-1^ 1/2 1/2 



dt. 



(2.23) 



where $(t) = $i + 1($2 - 4'i) with t e [0, 1]. 
For ah < a < 1/2, we have 

.1 .1 JV N 

/ /'*(.r'^^= / (E^?.+2^E^i 

•^0 i=l i=l 



AT 



/V^,2 , V^/ - - n2/, , E^Il </'l,i(<)^'2,i - 01,^)^^2 ( E^l ('/'2,j - (/>1,.)) 



2\ "-1/2 



,i=l i=l 



E.^l(<^2,i - 01,i)2 



dt 



< t + 



S»=l'/'l,»('/'2,. -'/'l,i) |2a-l 



E('?^2,i - (/-Li) 



2\"-l/2 ^ 1 Q-1/2 
Q,22a^*i-*2' 



which together with the fact that < p$(t) < 2(p$j + P p^.-^^^^.^) impUes that for all < a < 1 



n Jo 



/I I a-l\ 1/2 1/2 1/2 1/2 ^ , a/2 , x 1/2 1, 



1/2 



< lln"/2 



0,12/(6-Q),nl|Pr lloa2/(6-a),0 + ||P*i-$2 llo,3,0 WpTW 



.1/2 



1/2 I 



1/2 I 



< (11*1 -'i>2||?^o + 11*1 -*2|llo)ll*lli,o|ir||i. 

Similar arguments to that in (I2.2ip yield that 



(2.24) 



N 



E I^('?^1J'0J''^1:»'^) - -D((/'2jV'j>2,iW)| 



J=l 

AT 



N 



< E I^(01j'/'J ~ 02,i^/'j,(/'l,iW)| + E \^i^2,j'4'j,4'l,iV - 4'2,iV)\ 
j=l J = l 

AT AT 

< E ll'^^i'J ~ '?^2j||i,a||V'j||i,nlMli,o + E ll'^i.' ~ 02,»||i,n||V'jlli,n||w||i,n 
j=i j=i 

< ||*i-*2||i,a||*l!i,n||f||i,n \fi'en,yveH^{n). (2.25) 

Therefore, taking (EU), (EH]), (P?^ . (I^?^ . ([^TM)) and t^TI^i into account, we get 

iiPi - ^;)((n.*))> (2?.r)) < (llyi - y2\\Y + hi - 2;2||?^)||(p,*)||y||(x,r)||y v (p, (2^,r) e r, 

which implies p. 161) and completes the proof. 

□ 
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3 Finite dimensional approximations 

For the sake of generality, we will not concentrate on any specific approximation, rather we shall 
study approximations in a class of finite dimensional subspaces Sn X (n — 1,2, ■■■ ) that satisfy 

lim inf IIV^-^lli^o =0 V(/)eX, (3.1) 

where X is some Banach space containing the eigenfunctions of (j2.1[) . say, -ffo(il) or iJ^(il). 

Assumptions (j3.ip is apparently very mild and satisfied by several typical finite dimensional 
subspaces used in practice, for instance, spaces spanned by plane wave bases [7], spaces spanned by 
wavelets [51ll3|. and piecewise polynomial finite element spaces |llj . As a result, we may investigate 
all these kinds of finite dimensional approximation approaches in computational either physics or 
quantum chemistry in a unified framework. For convenience, here and hereafter we consider the 
case of X = i?o(Sl) only. 

We see that finite dimensional subspaces 

Hn = Sj^ c Ti 

satisfying 

lim inf |l*-«'|lin = V $ G 7^. (3.2) 

We shall study the numerical analysis of the following minimization problem 

inf{S($„):$„eH„nQ}. (3.3) 



The existence of a minimizer of p.3|) can be obtained by similar arguments to that in the proof of 
Thcorcm l4.1l (c.f.. also, [6l[8]). However, the uniqueness is unknown even up to a unitary transform. 
Therefore we define the set of finite dimensional ground state solutions: 



Gn - <^ e -«„ n Q : S($„) = min S(*) 
Given 71 > 1, any minimizer $„ = ((/)i,n, 02, «, • • ■ , 4'N,n) of ()3.3|) solves 

N 

(^$,>i,„,i;) = ['^Xij,n4'j.n,v) \/ V € Sn, i^l,2,---,N, 

(3.4) 



< 



In 

with the Lagrange multiplier 

N 



= {Kj.n)t j = l ^ { / 4>J,nA^„4>hn ] ■ (3.5) 
\Jn / ij=l 

Define the set of finite dimensional ground state eigenpairs 

e„ = {(A„, $„) e M^^^ X {-Hn n Q) : $„ e e„ and (A„, $„) solves ^} . 

Proposition [221 and p.Sp imply that the finite dimensional approximations are uniformly bounded 

sup (||*n||i,o + |A„|) <C (3.6) 
(An,*Ti)ee,i,n>i 
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for some constant C. 

We then address the Galerkin discretization of (|2.9I) . Let 

and -F„ : Yn — s- be an approximation of F defined by 

(F„((A„,$„)),(i?n,r„)) = (^^((A„,$„)),(x„,r„)) V (A„,$„),(x„,r„) e r„. 

Then discrete problem p.4l) can be rewritten as 

F„((A„,$„)) = Oei;*. (3.7) 
We also denote the derivative of F„ at (A„, <&„) e y„ by F'^ : F„ Y* as follows: 

AT 

{K,iA,^,9;^)iii^n,^n)),iln,'rn)) = (A„ , $„ )^'„ , r„) - ^ (^ij,„0j,n , 7i,n) 

iV 

Given (A, $) G 5^^^ x Q, we define 

X<j,=5^^^x(5$©r$)cr 
with the induced norm ||(|)J, ^) || = + for each {p,'^) G and 

a:$,„ = 5^^^ X (H„ n (5$ © r*)). 

We assume here and hereafter that yo = (Aq, $0) is a solution of p.5l) satisfying (|2.13p . where 
Ao = {'^0,13)^^=1 and $0 = (</'o,i' '^'0,2, • • • , 4'o.n)- We shah derive the existence of a unique local 
discrete solution y„ G -^£^$0,11 of ()3.4p in the neighborhood of yo- 

Lemma 3.1. Fy^ : X^^, — > Xl,^ is an isomorphism. 

Proof. It is sufficient to prove that equation 

= (3.8) 

is uniquely solvable in AT^^ for every (ir|,.g) G X^^. To this end we define the following bilinear 
forms a$, : 7^ X -H ^ M and 6*,,, c^^ : H x ]R^><^ M by 

a*„(vi/,r) = (/::,„(Ao,$o)*,r), 

i,J = l 

c<i>o(*!X) = ^ Xtj{i(l^o,i,ilJ]) + {(t)o,j,il'i))- 
Using (f2TT|) . we may rewrite ([3^ as follows: find |jj G 5^^^ and * G © Tio such that 

«$(,(*, r) -6$„(r,|jj) 

< 

CcE.o(*)X) 



(.9,r) vrG5$„©r$o, 

EV7 ^ cNxN (3-9) 
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For any given ^ £ gNxN ^ choose ^' = ^ol, and thus 

N 



c*„(vE',x) = 2 ^ , (3.10) 



where $q <i>o = jNxn jg used. Note that a simple calculation leads to 



N 

mi,n = ll^oxlli.n < iY. |x.,f )'/'ll<i>o||i,n. (3.11) 



By taking into account (|2.3p . p.lOp and p. lip , we obtain 



C$0 (^ 

XG5»x«^V5.^ ||*||l,n(E5=l IX^.I')'^' 



sup — ' , ^TTT^ ^ (3-12) 



where Kc > is independent of ^. Hence, there exists a unique solution 4*5 G iS$n such that 

N 



Therefore p.9p is equivalent to: find G 7io such that 

a$„(*o,r) = (g,r)-a$„(*s,r) vrGr$o. (3.13) 

The unique solvability of p.l3p is a direct consequence of p.l3p . 
Using similar arguments to that from p.lOp to p.l2p . we get 

mt sup > Kb, 

Xe5~x«^g5^^ ll*l|l,n(E5=llX..f)^/' " 

where > is independent of ^. This implies that equation 

6$„ (r, p) = a$„ (*o + *s, r) - (5, r) v r g 

has a unique solution £ S^^^ . 

We have proved that for any (111,5) £ P-9P - there exists a unique solution (jjs, ^'o + ^'s)- 

This indicates that Fy^ is an isomorphism from to and completes the proof. □ 

Note that Fy^ : X$q — ?> X^^ being an isomorphism is equivalent to the following inf-sup condition 

inf sup „ i;^^"^''^'^ -/3>0 (3.14) 

VieX^g y^f^X^g | | J/ 1 | | X^. j | ^ 2 | | X^. Q 

with the constant satisfying — \\Fy~-^\\. 
For any $ G Q, we define 

Q* = {*gQ: ||*-$||o.j2- min ||W-$||o.a}. 
In our analysis, we need the following lemma, whose proof is referred to [6]. 
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Lemma 3.2. // $ G Q, then 'J e Q* can be represented by 

= $ + 5(T^)$ + H^, 

where W e %> and S{W) G 5^^^ satisfying 

\S{W)\ = |(/^x^ ~ W^Wy/' - /^x^l < <\\^- <^>\\l^. (3.15) 

Before giving a discrete counterpart with Lemina l3.11 we also need to introduce two projections. 
First, we define the projection n„ : Q — > Hn H Q such that 

||n„*-*||i.fi = min V$e(Q. 

*e-H„nQ 

To project further into X$^„, we then define n„ : 5^^^ x Q ^ ^$.n by 

n„(A, $) = (A, (n„a>)c/) v (A, $) e 5^^^ x q, 

where 

;7 = arg min || (n„$)C/ - $||o,n- 
From Lemma 1321 we see that n„ : S'^^^ x Q ^ X$^„ is well-defined. 

Lemma 3.3. If Assumption A2 is satisfied, then there exists no > 1 such that F!^ n„j/o ■ "'^^o," ~^ 
n isomorphism for all n > uq. Moreover, there is a constant M > such that 

\\K,u^y,^\\<M Vn>no. 

Proof. We first prove that 

lim \\n,,y - y\\x, - V y = (A, $) g S^""^ x Q. (3.16) 

Using the fact that $ G Q and (n„<I>)?7 G Q*, we have 

1^7 - /| = \\iUn^)u - fi„$||o,o < ||(ri„$)?7 - $||o,j7 + ||n„$ - $||o,o < l|n„$ - 

which implies 

||(fi„$)?7-$||i,fi < ||n„$-$||i,fi + ||(n„$)J7-n„$||i.o 

< ||fi„$-$||i,o + |c/-/| • l|n„$lli,o 

< ||fl„$-$||i,o. (3.17) 
Let = (01 , 02 I • • • I '/'tv) = arg min - we may estimate ||n„<i> - ^Hi.n as follows: 



N 



t IIQ™dlo,n 



Pi\\l,Q. 



N 



^ \\Qn(p"\\o,n 



i=l 
N 
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where Q„ is the Gram-Schmidt orthogonal operator: 



Note that 



^ WWnfPj Ho o 



3 = 



we conclude 



^ /I , IIQn'^^i lll.S^Mi , ,„,, , ^ IIQn'/'jIll.n II , „ ,„|| 



||n„$-$||i.f2< inf ll^-^lli.n. (3.18) 



Using (I3.17p . (|3.18l) and the definition of n„, we arrive at 

mny-y\\x, < inf ||*-<i>||i,n, (3.19) 

which together with p.2p leads to p.l6p . 

We then show the invertibility of F,' : X<j,g^n ~^ „. We obtain from p.l4p that 

sup -r II II >li V yi e A$o,„. 

y2eXi„ \\yi\\x^^ 112/211X4,0 

Let P*° : 5<5(, n %>„ T~in H (<S$o fl 7$^) be a projection operator satisfying 

(V$i,V($2-nf"$2)) =0 V$i eH„n(5*„nr*o). 

Set 

sup ||^'-P*°^'l|o,n, 

*e54,Qnr*(,,||*||i,s-!<i 

we have (see, e.g., [34] ) 

||*-F*°*||o,o <??„I|*||i,o V*e5$„n7io with lim?7„ = 0. (3.20) 

n^oo 

Let f„ = (/, .^ye obtain from definition (piTj) and that 

{F^^yuPny2) - (F^„2/i,j/2)-(^^;„2/i,y2-P„y2) 

= (F^o2/i,y2) + ^(V<i>i, v($2 - P*"$2)) - (i^^,yi,y2 - P„y2) 



> {Fy^yi,y2) - C||yi||x*j|y2 - ^'n2/2||o. 



which implies that there exists fi such that for all n > ri, there holds 

sup T. — rr-n — \\ ^77 v G ^$0,") 



y2eJ>s:4.o,„ ll2/i|k«.Q I|y2||x 



*0 
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or equivalently 

mf sup — — ^ > — . 

yi6^*o."a2ex*„,„ 112/1 ||x^J|2;2|U4,„ 2 

Thus F^^y^ is an isomorphism from to „ satisfying 

IK,;' II < 2/3-' vn>n. 

Note that satisfies the following discrete Holder condition 

IK,,„ - <,n„,„ll < llyo - n„yoll^,„ + bo - n„2/ollL„- 

It follows from p.l6p that there exists no > n such that the inf-sup constant of F'^ iinyo uniformly 
away from zero for all n > uq. This completes the proof. □ 

Theorem 3.1. If Assumption A2 is satisfied, then there exist 5 > 0, ni > uq such that p.4p has 
a unique local solution j/„ = (A„, $„) £ X$o-n ^ Bs{yo) for all n > ni. 

Proof. The idea is to construct a contractive mapping whose fixed point is We rewrite p.7p as 

FniVn) - F.niU,,yo) - -F„iU„yo). 

Using (|2.15p . we have 

||K(n„yo)IUi^,„ - ll^^(n„2/o)|x,„,„-^^(2/o)|x,„,J|x;^,„ 

< ||F(n„yo)-F(yo)|ki^^ < ||2/o-n„yo||x,„. 

From Lemma [3.31 we may define the map TV : i3ii;(n„j/o) H X'S>o,n X<s>o,n by 

Jo 

when n > hq. 

We will show that A/" is a contraction from i3fl;(n„yo) H A'$(,,n into Bjf (n„yo) H A^^o," if i^ 
chosen sufficiently small and n is large enough. 

First, we prove that J\f maps i?fl(n„2/o) H ^*o," to Sfl'(n„?/o) H X,pg.n for sufficiently small i?. 
Note that F^ ii„yo isomorphism on X$,-|.„ if n is sufficiently large. For each x G i?7j(n„?/o), we 
have N'{x) — n„?/o G -''^^o," ^^^^ 

||AA(x)-n„2/o||x,„ 

< M(||i^„(n„yo)IUi^,„+i?^ \\FU.yo+ti.-n^yo) - KM.yoWdt) 

< CM{\\n„yo-yo\\x^„+R{R" + R^))- 

Since CM(l|n„?/o — yollx^,;, + R^^" + R'^) can be estimated by R when R is sufficiently small and n 
is sufficiently large, we have that Af{x) G i3i^(n„yo)- It is clear that R can be chosen independently 
of n. 

Next, we show that A/" is a contraction on -B_R(n„j/o) H If xi,X2 G Si<;(n„yo) H 

then 

K,n.yo(^(^^)-^(^2)) = (^1 /' (^^;n„,o 

Jo 
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Thus, l|A/'(xi) — N'{x2)\\x^g can be estimated as 

||AA(xi)-AA(x2)||x.„ 

< M\\x2-Xi\\x^„ \\K,nr,yo ^ K.xi+t(x2-x^)\\dt 

Jo 

< CM{R" + R^)\\xi-X2\\x^^. 

We obtain for sufficiently small R that CM{R°' + R^) < 1 and hence A/" is a contraction on 
BniUnyo)- 

We are now able to use Banach's Fixed Point Theorem to obtain the existence and uniqueness 
of a fixed point j/„ of map Af : -B_R(n„T/o) H -^^$0," ~^ -Bii:(n„yo) H X^^^n, which is the solution of 
FniUn) ~ 0. This completes the proof. □ 



4 Numerical analysis 

In this section, we shall prove the convergence of finite dimensional approximations and derive 
various error estimates under different assumptions. 



4.1 Convergence 

The purpose of this subsection is to prove the convergence of the numerical ground state solutions, 
for which we need to introduce the following distances between two sets. We define the distance 
between two subsets A, B C Y hy 

V{A,B)^ sup inf (|A-p| + ||$-«'||i,a) 

(A,*)eA (iJ>*)eB 

and the distance between two sets M, N C H hy 

dH{M,N)^ sup inf 

Theorem 4.1. There hold 



lim D(e„,e) = 0, (4.1) 



lim = min (4.2) 



where En — i?($n) for any <!>„ g Q„ 



Proof. Let (A„, <&„) € 0„ for n = 1, 2, • ■ • . Given any subsequence of with 1 < ni < 
712 < • • • < Tifc < • • • , we obtain from the Banach-Alaoglu Theorem and (13. 6p that there exist $ G H 
and a weakly convergent subsequence {^uk- } <^ {'^'nt} such that 

$„,^ ^ $ in H. (4.3) 

Next we shall prove ^ € Q and 

lim II $-$11 1,^,-0, (4.4) 

lim £:($„, ) = min £:(«'). (4.5) 
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From (|4.3p and Proposition 12. 3[ we have 

lim inf . ) > . (4.6) 



J-S-OO 



Note that (|3.2I) implies that {^nt } is a minimizing sequence for E{'i') and the ReUich theorem 
shows that 



n Jit 
Therefore $ G Q is a minimizer of which together with (|4.6p leads to 

lim ) = Em = mini?(vl/). (4.7) 

This further implies (|4.5p and $ £ 

Since -ffg (f^) is compactly imbedded into L'^iyt) for p £ [2, 6), we have that 0i,„^. (f>i strongly 
in L^'(ri) as j — > oo for i = 1,2, •• • ,N. This indicates that {p<s>„ } converges to p$ strongly in 
L'^(ri) for q G [1,3), from which wc obtain that 

lim / Vioc{x){p,t,„ {x) ~ p<s,{x))dx ^ 0, 
lim / (f (p$ ) — £{p^{x)))dx — 0, 

and 

lim ) = D(p$,p$). (4.8) 

Consequently, we can get from (|4.7p to (|4.8I) that each term of £^(-) converges and in particular 

N N 



— >oo — ^ J ' — ^ 



7— >oo 

1=1 



/AT \ 1/2 

Using (|4.3p and the fact that ?^ is a Hilbert space under norm j ||V(/)i||o q j , we obtain (|4.4p . 
If (A, $) solves (USD, then 

hm |A - A„ I = 

is a direct consequence of (|2.7p . (|3.5p and (|4.4p . Hence we arrive at (|4.ip . This completes the 
proof. □ 

Remark 4.1. Theorem \4-l\ states that all the limit points of finite dimensional approximations 
are ground state solutions. We note that \2S^ gave the convergence of ground state energy ap- 
proximations only while we provide further convergence of approximations of both eigenvalues and 
eigenf unctions. 
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4.2 Error estimates for the energy approximation 

We shall derive the quadratic convergence rate of ground state energy approximations, which is a 
generalization and improvement of [SJ [55] • 



Theorem 4.2. Let E he the ground state energy of (|2.4p and En he the ground state energy of 
p.3p . namely, E = E{<^) for all ^ £ Q and En = for all £ If Assumption Al holds, 

then 

\E-En\<dl^{g,nn). (4.9) 

Proof. We see from the definition of ground state energies E and En that 

Q<En-E < E{^) -E y^eUnnQ. 

Following [HI [57], if Assumption Al holds, we obtain from the Taylor expansion that for any 
* G Q, there holds 

- E{<i>) = - $) + i(£;"(0(* (4.10) 

where ^ = <i> + S{'i> — $) with S E [0, 1]. Since $ is a ground state solution, we get from (|2.5p that 

(£:'($), * - $) = 2($A, * - $) = 2($;7f7^Af7, *c/ - m), 

where the orthogonal transform U diagonalizes the Lagrange multiplier A by 

C/^AC/ = diag{Ai,-- - ,~Xn}. 
Denote I* = $C/ and ^ = , we have 

N ^ N 



i=i ■'^ .=1 
< ll^-*ll?o = ll$-*ll?o- (4.11) 



It is observed by a simple calculation that 

N N ^ r 

(£;"(C)f , r) = 2 ^(A^V^, 7^) + 4 ^ i?(C^V'», e,7j) + 4 ^ / £"(p«)e.V'.07, 

and hence 

(i?"(0(* - a>), * - $) < - $112 (4.12) 

where the hidden constant depends on the H-norm of vf. 

Taking (|4.10p . (14. lip and (|4.12p into account, we have proved that for $ e there holds 

E{-^)-E{<P)<\\<P-^\\ln V*eH„nQ, 

which together with the definition of n„ and (|3.18p implies that n„<i> G Hn H Q and 
0<En-E< E{fln<i>) - S($) < ||n„$ - nln < dUG^nn), 



where the hidden constant, by using p.6p . is only dependent on the problem. This completes the 
proof. □ 
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4.3 Error estimates for ground state solutions 

In this subsection, we shall derive the a priori error estimates for finite dimensional approximations 
of Kohn-Sham equations under Assumptions A2 and A3. Note that yo = (Aq, $o) is a solution of 
(|23|) satisfying ((2l3t . 

We define bilinear form a'(<i>o; •, •) by 

a'($o; r) = (Ao, $o)*, r) v r e 

Obviously, a'($o; •) is continuous on "H x Ti. 

Now we shall introduce the following adjoint problem: for / e (L^(f2))^, find ^/ £ 7*0 such 
that 

a'($o;*/,r) = (/,r) vrer*„. (4.13) 

Since C'^^{Aq, $0) is an isomorphism, (I4.13|) has a unique solution and 

ll*/lli,n< ||/||o,n. (4.14) 
Let K : {{L'^{n))'^, (•, •)) ^ (T$o, (V-, V-)) be the operator satisfying 

{VKw,vv) ^ {w,v) y w eiL^in))'^,y V (4.15) 

Then K is compact. Set 

Pn= sup inf ||V((/:'^„(Ao,$o))-^X/-*)||o,a, 

/e(L2(n))«,||/||o,n<l*e^" 

we then have the following estimate (see, e.g., [4]) 

||V((/:^„(Ao,$o))-iif/-P;(/::,„(Ao,$o))-iif/)||o,o<P„||/||o.o yfeiL'mf (4.16) 
with 

lim p„ = 0, 

n— >C30 

where : Tisig 7io H H„ is the projection operator satisfying 

(V($i - P>i), V$2) = V $2 e r*„ n H„. 

Theorem 4.3. // Assumptions A2 anrf A3 are satisfied, then there exists 5 > such that for 
sufficiently large n, p.4p has a unique local solution (A„, $„) G Ar$o^„ n Bs{yo) satisfying 

||*o-$«||i,f2<d«(e,H„) (4.17) 

and 

11*0 - $«||o,f2 + |Ao - A„| < p„||$o - *n||i,a (4.18) 

with pn as n oo. 

Proof. We obtain from Theorem 13. II that there exists S > such that for sufficiently large n, (13. 4p 
has a unique local solution y„ = (A„, $„) e Ar$(,^„ n Bs{yo). Hence, we have 

FniVn) - -Fn(n„yo) = --Fn(n„yo), 
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which leads to 

K.n„y„iyn - n„yo) = -i^„(n„yo) - iVn - n„2/o) / (-F'«,n„yo+t(y„-n„yo) - K,n^yo)dt- 

Jo 

Using the similar arguments in the proof of Theorem I3.1[ we obtain from Lemma 13.31 that for 
sufficiently large n 

-n„?;olU^o - llyo -n„?/o||x^„ + hn - 'nnyo\\x^„i\\yn - 'nnVoWx^^ + bn -n„yo||x.^J, 

which together with p.l6p and the fact that ?/„ G Bgd/o) implies that for sufficiently large n 

||2/„-n„2/o||x^„ < bo-n„yolU^„- (4.19) 
Using (IXTg|) and we conclude 

hn - yo\\x^,, < hn - n„2/olk*o + \\yo - ^nVoWx^g < inf 11^- - ^olIl.O, 

which implies (|4.17p . 

Since there exists Si G [0, 1] such that 

Jn 

where ^ = (Ci, 6, • • • , ^n) with = (5i0i^„ + (1 ^ <5i)0o,i, using Assumption A2 we get 

^ llo,3/a,o||</'i,n " 0O,i || 0,f2 1 1 (/"j-.n || 0,6/(3-2a),O + 1 1 Pf 1 1 0,3,^ || " 0O,i 1 1 0,0 || 0j,n || 0,6,n 

< \\4>i.n - 4>o,ih.n, 
from which we have 

((^'(P*,J - ^'(P*o))'/'»,"'<^j:n) 
= (^^'(P*„ n - £'{P<i'a)4>0,iA],n) + ' (p*o ) ('/'O.i " 4>i,n), 4'j,n) 



Note that 



^ n 

= (A$o((?!>i,n - 00,*), - <l>0,j) + y2 >^0,ik<t>0.k{4'j,n " <t>0.]) 

■1^ k=l 

+ / Aojfe</'O,fe(0j,n - 0O,j) + / (^'(P*„) ~ ^'(P*o))'/'»,n<^J.n 

+D{(t>i,n4'j,n, P<S>„ -P$o)- 

Hence we conclude that 

|A„ - Aol < ||*„ - $o||?,o + - 3'o||o,n. (4.20) 
By Lemma [321 we decompose $„ as 

= $0 + 5(W^)4'o + Vr, (4.21) 
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where e 7io and S{W) e 5^^^^ satisfymg 

\S{W)\ < \\W\\l,,<\\<i>o~<^n\\ln- (4.22) 
Setting 5' — and applying the duahty problem of (I4.13p . we obtain 

\\^n-^o\\l,n = («>„ - $o,*« - *o) 

= ($„-$o,5(VK)$o) + ($«-$o,VF) 

= ($„-$o,5(VK)$o)+a'(*o;*,VK), 

which together with (|4.2ip leads to 

- *o|lafi - - *o, 5(VF)$o) - a'($o; 5(VF)$o) + a'(*o; - *o) 

= ($„ - $0, 5(W^)$o) - a'(*o; 5(M^)$o) + a'($o; * - K'^, - *o) 
+a'($o;-P;«',*„-$o). 
Note that from (|2.5p and p.4p . we have 



while the fact that G Ti,, yields 

Jn Jn Jn 

we then come to 

ll*n - $o|lo.n - - $o,'5(VF)$o) - a'($o; ^',5(VF)$o) + a'($o; * - ^l*, - $o) 

-i(i?'(ci>„)(P^5') - i?'($o)(P^M') - i?"($o)(P>, - $o)) 

Using the Taylor expansion, we have that there exists S G [0, 1] such that 

i?'($„)(F^vl,) _ E'i<S>o){K^) - E"i<S>o){P:^^,<S>^ - $0) 
= E"iOiP:,^i!, - $0) - i?"($o)(P^*, <!>„ - <i>o) 

< (ll<i>n-'i>o||?,n + ll*«-$o||?,o)ll<fn-*o||^,o, (4.23) 

where ^ = $0 + <^(*i')i — ^0) and the last inequality is obtained by the similar arguments in the proof 
of (|2.22p or Lemma 4.5 in 6 when Fi = $„ - $0, ^2 = ^71 - ^0 and = P,' and using the fact 

ll^'n^lli.n < ||*||i,n < ||<i>n-<i>o||o.n. 
Taking (|4T4l) . (|4?T6l) . (|422l) and (|4?23l) into account, we obtain that 

||$« - *o||o,n < ||$„ - $o||g.o + PnW^n ' *o||i,n + ||$„ - $o||"nll*n - *o||o,n 
+ |A„-Ao|(||$„-$o||o,n+PrO: 
which together with (I4.20p and Theorem 14. II produces 

- $o||o,n < Pn||$« - $o||i,n 
when n ^ 1. This completes the proof. □ 
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Remark 4.2. Theorem shows that under certain assumptions every ground state solution can 
be approximated with some convergent rate by finite dimensional solutions. We see that ^61 pro- 
vided numerical analysis of plane wave approximations only while our results apply to general finite 
dimensional discretizations and the analysis is systematic and carried out under very mild assump- 
tions. 

Remark 4.3. If in addition, Vioc G H^{^), Cj S H^{n) {j = 1,2,--- ,M) and £ € Ci([0,oo)) n 
C'^((0, oo)), then for sufficiently large n, estimates (I4.17P and (I4.18P are also satisfied with pn ~^ 
as n — i> oo. Here 

p„= sup inf ||V((£'<,,,(Ao,<i>o))-ii^/-vI/)||o,o 

/6«,ll/l|i,n<l*«^^" 

andK: (H,(V-,V-)) ^ (r$o , (V-, V-)) .satisfying BIB . 

Remark 4.4. Let yo = (Aq, $o) be the ground state solution of (|2.5p satisfying (|2.13p . We assume 
that is a convex bounded domain and Sn is replaced by the standard finite element space SQ''''{fl) 
of piecewise polynomials of degree fc (fc = 1, 2) of HQ{fl) over a shape-regular mesh with size h. Let 
{Afi^k,^h,k) G ^io,h be the ground state solution of (j3.4p and Assumption A2 hold. Then 

|Ao - A,,,i| + ||$o - */i,illo,n + h\\<pQ - < h"^ 

when /i < 1. If in addition, Vioc G H^{n), Q £ H^{n) {j = 1,2,-- - ,Af) and £ e C\[0,oo))n 
C3((0,cx))), then 

|Ao - Ah,2\ + hW^o ~ ^h,2\\o,n + /i^||*o - $h,2||i,o < 

when h <^ I. 

5 Numerical examples 

In this section, we will report several numerical examples that support our theory. These numerical 
experiments were carried out on LSSC3 cluster in the State Key Laboratory of Scientific and 
Engineering Computing, Chinese Academy of Sciences. Our code is based on the PHG finite 
element toolbox developed in the State Key Laboratory of Scientific and Engineering Computing, 
Chinese Academy of Sciences. 

In these examples, we solved Kohn-Sham equation (|2.5I) . We chose our computational domain 
n as [—10.0, 10.0]'^. In computation, we used the norm-conserving pseudopotential ^SOj obtained by 
fhi98PP software and applied the local density approximation (LDA) for the exchange-correction 
potential. We applied the standard linear and quadratic finite element discretizations over uniform 
tetrahedral triangulations. The finite dimensional nonlinear eigenvalue problems were then solved 
by self consistent field iterations. In each iteration, the Kohn-Sham Hamiltonian is constructed 
from a trial electron density, the electron density is then obtained from the low-lying eigenfunctions 
of the discretized Hamiltonian, the resulting electron density and the trial electron density are 
then mixed and form a new trial electron density. The loop continues until self-consistency of the 
electron density is reached. 

We present numerical results for N2, C2II4 and Sill4 molecules. Since analytical solutions are 
not available, we use the numerical solutions on a very fine grid for references to calculate the 
approximation errors. 

Let us first come to the ground state total energy approximations. The errors of total energy 
of N2, C2II4 and Sill4 are presented in Figures [5TT1 15.21 and l5.3i respectively. We can see that 
convergence rates for linear and quadratic finite elements are h^ and h^ respectively, which agrees 



22 



well with the results predicted by Theorem 14.21 We then present the approximation errors of the 
first two eigenvalues for these three molecules, see Figures 15.51 and 15.61 We may see that these 
results coincide well with our theory (see, e.g.. Remark 14.41) . too. 





10'"*^ 10""^ 10"°^ 10° 10°^ lO'*''^ lO'"" 10"°^ 10° 10°^ 

mesh size h mesh size h 

(a) Linear finite elements (b) Quadratic finite elements 

Figure 5.2: C2H4: errors of the ground state total energy 



6 Concluding remarks 

We have analyzed finite dimensional approximations of Kohn-Sham models. We have proved the 
convergence and shown the optimal a priori error estimates of finite dimensional approximations. 

As we see, the ground state solutions oscillate near the nuclei |14[ ITT] . It is natural to apply 
adaptive finite element discretizations to carry out the electronic structure calculations. Indeed, it 
is our on-going work to study the convergence and complexity of adaptive finite element methods 
that will be addressed elsewhere. 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 5.3: SiHi: errors of the ground state total energy 
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(a) Linear finite elements (b) Quadratic finite elements 

Figure 5.4: N2: errors of the first and second eigenvalues 
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Figure 5.5: C2H4: errors of the first and second eigenvalues 
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